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Abstract 

We consider primordial spectra with simple power behaviours and show that in the Navier- 
Stokes and magnetohydrodynamics equations without forcing, there exists systems in three 
dimensions with a subsequent inverse cascade, transferring energy from small to large spatial 
scales. This can have consequences in astrophysics for the evolution of density fluctuations, for 
primordial magnetic fields, and for the effect of diffusion. In general, if the initial spectrum is 
k a , then in the "inertia!" range, for a > —3 there is an inverse cascade, whereas for a < —3 
there is a forward cascade. 
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Recently there has been a considerable discussion of the development of primordial magnetic 
fields subsequent to the creation of these fields (see e.g. refs. JTJ, 0, and ||, where further 
references can be found). In the early universe the Reynold number Re is often quite large. Since 
Re is of the order the ratio of the nonlinear to the linear terms, this suggest that the nonlinearities 
could be quite important. This was indeed found in ref. |2| and ||, where an "inverse cascade" 
was observed. In this note we point out that it is an exact consequence of the Navier-Stokes 
and magnetohydrodynamics equations without external forces that there exists an inverse cascade, 
moving the field to larger scales, provided the initial energy is distributed at relatively small scales. 
This result can be expressed as a scaling relation for the relevant energy. 

In astrophysics one often considers "primordial" spectra behaving like some power k a . The 
question we address in the following is how such spectra evolve as consequence of the non-linear 
equations of motion. We shall find that if the such a system is "left alone" (like in many astrophysics 
applications), its energy is moved to larger and larger scales, in contrast to a "forced" system 
(relevant for earthbound hydrodynamics) in three dimensions, where energy is moved from larger 
to smaller scales. Hence there is a fundamental difference between systems which are forced or left 
alone. 

In this note we shall consider hydrodynamics or magnetohydrodynamics (MHD). In hydrody- 
namics one often considers systems which are under the influence of forces, e.g. at large distances 
("stirring"). In these systems there is a forward cascade in three dimensions, which means that 
energy is transferred from large scales to smaller scales (from "order to chaos"). As mentioned 
before, cases without forcing may also be of interest. In particular, in astrophysics one may of- 
ten encounter systems without any essential "stirring" , where the initial "primordial" spectrum is 
given. It is then quite dangerous to apply results from forced systems. Here we shall show that for 
a certain initial spectrum there exist an inverse cascade in hydrodynamics (from "chaos to order"), 
whereby energy is transferred from smaller scales to larger scales, as a consequence of the exact 
Navier-Stokes or magnetohydrodynamics (MHD) equations. 

The most rigorous of our results is that if the primordial energy spectrum is given by k at time 
0, then at later times it will evolve as 

k ^J(k 2 t), (1) 

where tp is a scaling function (provided the viscosity coefficient is constant in time) with ^(0) = 1. 
For MHD an equally rigorous result can be derived for the magnetic energy. These results take 
into account diffusion. In the "inertial" range, where diffusion can be ignored, a much more general 
result, relevant for an initial spectrum k a , can be derived. In order to have a finite energy, we must 
introduce an ultraviolet cutoff to be specified more precisely later, so that k < K. 

In a more special context the scaling behavior (|l]) of the energy density has been considered 
before, first by Heisenberg j| in his effective diffusion model, and then by Parisi || in a continuous 
version of the cascade (GOY) model ||. A generalization of the scaling in (|l|) was then found 
in an investigation of the inverse cascade in the continuous GOY version of three dimensional 
relativistic MHD by Brandenburg et al. ||. Later the scaling (|l]) has also been generalized to 
turbulent mixtures In the following the result (|l]) and its generalization to MHD are derived 
from the exact Navier-Stokes and MHD equations, showing that Eq. (|l]) is much more general than 
previously thought []. 

1 This result can be used to check if model calculations are in agreement with the Navier-Stokes equation: If 
forcing is taken to be absent, the scaling (|l|) should be valid. From this point of view the cascade (GOY) models are 
excellent, since they have the scaling as an exact solution in the force-free case, as shown in Refs. p), H, and M. 
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We start from the expression for the energy in k = |k| -space, 



E(k, t,L,K) = —- d s xd 3 y e^^) < v(x, t) v(y, t) >, (2) 



(2k) 

with the total energy given by 

dk E(k, t, L, K) = / 

l2n/L J2tt/K 



[ K dk E(k,t,L,K) = [ L d 3 x - < v(x,t) 2 >, (3) 

hitlL J2w/K 2 



where L indicates the insertion of an infrared cutoff to be specified precisely later. 

Our initial spectra are taken to be simple powers. Thus, for powers larger than or equal -1, we 
need the ultraviolet cutoff K mentioned before. If this is translated to the behaviour of the velocity 
correlator in the integrand on the right hand side of Eq. (|2|), we need a short distance cutoff 1/K 
in x-space at time t = 0. This is because the primordial spectrum is generated by some physical 
mechanism, and hence the initial energy must be finite. Therefore K is determined by the physical 
mechanism which generates the initial spectrum. For simplicity we assume that the cutoff can be 
represented by a factor^. 

F(k/K), where F(x) ~ 1 for x < 1, 

and F(x) « for x > 1, (4) 

for t = in the primordial energy E. When t increases, the initial cutoff is no longer needed 
because of diffusion, from which we expect the energy to decrease like exp(— 2vk 2 t) for large values 
of k. 

Next we want to use the well known self-similarity property of the non-relativistic Navier-Stokes 
or MHD-equations, 

x -» Ix, t -» l l ~ h t, v -» l h v, v -» l 1+h u, B -» l h B, 77 -» Z 1+/l ry, (5) 

where i/ is the kinetic and r; is the Ohmic diffusion. Using the substitutions x = lx' and y = Zy', 
we obtain from eqs. (Q) and (||) 

E(k/l,l l - h t,Ll,K/l) = l 4 ^^ dVdV e' k(x '" y,) < vf/x',/ 1 -^) vf/y',/ 1 -^) > 

(27T) J J2tt/X 

= l i+2h E{k,t,L,K). (6) 

Here 2£ deviates from E 1 in the sense that it is the energy computed from a theory where v is 
rescaled according to (|5|). The cutoff K, introduced by (Q), has the property from (||) that it only 
depend on k/K, and hence is the same factor on both sides of Eq. (||). For t > viscosity provides 
the cutoff by the factor exp(— 2vk 2 t) in the energy. However, the quantity vk 2 t is invariant under 
the self-similarity transformation (|5|), and hence the diffusive cutoff is also the same on both sides 
of Eq. (g). Thus, for k < K we can leave out the independence on both sides of Eq. (||), and in 
the following we therefore leave out the argument ll K v in the energy, and we tacitly assume that 
k < K in the following. 

Let us define as usual the "inertial range" as the interval of /c-values where the dynamics is 
independent of the scale of dissipation, so that in this range diffusivity can be taken to be zero. 
We assume that such a range of k- values actually exists. If one of the conditions 

2 For example, if F(x) = e~ x , then this can be achieved in x-space by replacing some factors of the type l/(x — y) 2 
by combinations of factors the type l/(jx — y| ± i/K). 
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(i) h = -l 



(ii) h 7^ — 1, but we are in the "inertial" range, 

is satisfied, we have E(k, t, L) = E(k, t, L). This is because if (i) is satisfied, diffusion is left invariant 
in the self-similarity (||), and if (ii) is satisfied, we can take z^=0. In the following it is assumed 
that one of these conditions are satisfied, and we therefore take E = E. Eq. (|6|) then becomes 
a functional relation for the energy, which we shall solve in the following. Introducing the energy 
density £ by E = £V, where V is the normalization volume, defined with the infrared cutoff as 

V = J^d 3 x, (7) 

we have 

£{k/l, l^H, IL) = l 1+2h £(k, t, L). (8) 

If the infrared cutoff is implemented in x-space by an exponential exp(— |x— y|/L) in the integrals in 
Eqs. (Q) and (^), it will appear in /c-space by the substitution k 2 — > k 2 + 1/L 2 for some factors of k 
in the spectrum £. This is because the exponential cutoff, after performing the angular integration, 
is equivalent to the substitution |k| — > |k| ± i/L in the integrals in (Q) and (|6|). For the simple 
primordial power spectra used in this note, it therefore follows that as far as £ is concerned, the 
limit L — > oo can be performed without encountering any singularities except possibly in the point 
k = 0, i.e. 

£(k,t) = lim £{k,t,L) (9) 

exists. Hence the cutoff dependence in Eq. (||) can be ignored in the following. 

The exception to Eq. (|9|) is that if £ is formally singular in k = 0. In order to have the 
physically required finite energy, k should be substituted by y/k' 2 + 1/L' 2 . The value of L should 
then be determined by the physics which generates the primordial energy. 

If we introduce the function ip 

£{k,t) = k- 1 ~ 2h i){k,t), (10) 

we see that 

^{k/l,l l - h t) = ^{k,t). (11) 



Differentiating (11) with respect to / and putting 1=1 afterwards, we get 



with the solution ip = /((l — ft) log A; + logi), where / is an arbitrary function. Thus, using eq. 
([l0|) we have the scaling relation 

£{k,t) = k- 1 - 2h -ifjik 1 -^). (13) 

As is clear from the preceeding discussion, this relation is valid if either h = —1, or, more generally, 
if we are in the inertial range where diffusion can be ignored. 

The general interpretation of ((l^) is that if initially the "primordial" spectrum is given by 

£{k,0) = k a , (14) 
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then at later times the spectrum in the inertial range is given by 

£{k,t) = k Q V(A; {3+Q)/2 t). (15) 

Depending on the value of a we may need an infrared cutoff. If h = —1, we see from the self- 
similarity (||) that if the initial spectrum has a=l {h = —1), then the spectrum is always given by 



eq. (15), because the diffusion coefficient is the same in the scaled and unsealed equations, i.e. 

£(k,t) =kiP(k 2 t) (16) 

holds as a consequence of the Navier-Stokes equations for all values of k. It should be noticed that 
the reason this type of scaling works also when diffusion is included, is that diffusion scales exactly 
like in Eq. (|i~6|). The new thing is that the non-linear effects scale the same way, if the initial 
spectrum is linear in k. Thus we see an inverse cascade, instead of a cascade^]. 

In the case of MHD, Eq. (|l6|) is still valid, and is to be supplemented by the relation 

£ B (k,t) = k (j)(k 2 t), (17) 

where £b is the magnetic energy density, 

S B {k,t) = J d s xd 3 y e ik < x -y) < B(x,i)B(y,t) >, (18) 

and 4> is another scaling function. We also used that B scales the same way as v, and the Ohmic 
diffusion scales like v. It is now quite clear that the spectrum has an inverse cascade: The spectra 
start at t=0 with the linear form £, £b oc k, and then the dynamics at later times is governed by 
the scaling variable k 2 t, which is invariant for k 2 oc 1/t, so at later times, for a given value of tp, k 
diminishes, and energy is transferred from smaller to larger scales. 

From Eq. ( p"5|) we see that in general in the inertial range there is an inverse cascade as long 
as a > —3, otherwise we have a direct cascade. This implies that if a (magneto-)hydrodynamical 
system is left to itself, the way in which it develops is governed by the initial distribution of the 
energy (for k < K): If much initial energy is concentrated at large (small) distances, the system 
will cascade towards small (large) distances. 

So far, we have considered the non-relativistic (magneto-)hydrodynamics equations. The general 
relativistic case for an expanding, flat universe has been discussed in ref. [g]. The result is that the 
equations have the same form as in a non-expanding universe, provided we scale the energy density 
by i? 4 and the magnetic field by R 2 , where R is the scale factor. Furthermore, the dynamical time 
now becomes conformal time t, defined by 

t = j dt'/R(t') oc R(t) oc Vi. (19) 

Here t is the Hubble time. The dynamical time t is thus slowed down relative to the Hubble time. 
The velocity is not scaled, and the coordinates entering differentiation are the comoving ones. 
Therefore, k is the comoving wave vector. The conserved energy is the energy multiplied by i? 4 . 
Therefore Eq. ( |l5j ) is changed to 

£(k, t)R(t) 4 = k a -0(A; {3+Q)/2 t~) = k a 7p(R(t)k {3+a ^ 2 ). (20) 



3 It goes without saying that if if) is dominated by diffusion, the result is not interesting. Therefore the Reynold 
number should be so large that non-linear terms are important. 
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Also, the magnetic energy is modified to 

S B {k, t)R(t) 4 = k a 4>{k^ +a ^ 2 t) = k a 4>(R{t)k (3+a ^ 2 ). (21) 

These expressions are valid if assumption (ii) is satisfied, i.e. in the inertial range where one can 
ignore diffusion. In most cases of interest in astrophysics, diffusion is not constant in time. For 
example, for Silk damping the effective diffusion behaves as v/R oc t 2 . In such a case we can 
only say that the inverse cascade competes with diffusion. Detailed calculations at high Reynold 
numbers are needed in order to see whether the inverse cascade wins over the diffusion (see Ref. S). 
It should be noticed that the results for the expanding universe are valid in comoving coordinates. 
The expansion of the physical coordinates, or the red shift of k, comes on top of the inverse cascade. 

The fact that a (magneto-) hydrodynamic system behaves in a way which is fundamentally 
different when it is forced or left alone, is perhaps not of much interest in the usual earthbound 
hydrodynamics, since systems without forcing die out. However, it can have consequences for 
astrophysics, because the various spectra, e.g. density fluctuations, observed now may have been 
subject to inverse cascades earlier. Thus, going backwards in time, the primordial scales may have 
been smaller than what is generally thought from data obtained now. This is especially true since 
in general the relevant Reynold numbers are quite large in cosmology, so the non-linear terms are 
expected to be important. Of course, in order to get quantitative estimates, it is necessary to do 
explicit calculations in order to determine the spectra, or equivalently, the scaling functions ip and 
<j). This was already done for the magnetic energy spectra by Brandenburg et al. in a shell model 
for relativistic MHD in a flat, expanding universe, see Refs. § and Q. Depending on the actual 
magnitude of the inverse cascade, it may, for example, imply that less inflation is needed. 

An inverse cascade has the effect of generating order from chaos. This was clearly seen in 
Ref. [||. To give an example, in the case where the initial spectrum is linear in k, this spectrum 
corresponds to initial short range velocity correlations in x-space. However, as time passes, the 
correlation gets a longer range, because k oc l/\/t, thereby creating more order. 

In the previous discussion we have ignored the effect of gravity (except for the inclusion of the 
scale factor). Since gravity is an attractive force, it will transfer energy from larger to smaller scales, 
and hence it will compete with the inverse cascade. In the non-relativistic case there is, however, 
one case where there is still scaling with gravity included, namely when h = —1/2. In this special 
case the gravitational force scales like the other terms in the Navier-Stokes equation. For h = —1/2 
one has the initial energy £(k,0) =const. In this case, at later times one has £(k,t) = i;(k 3 / 2 t) in 
the inertial range. 

Finally, inverse cascades have the general effect of diminishing diffusivity, because energy is 
transferred away from the small scale region, where diffusion is operative. Therefore, there will be 
more energy at larger scales than one might infer by ignoring the non-linear terms in the equations 
of motion. For MHD, the diminished effect of Silk damping Q was seen in Ref. |J. 

I thank Larry McLerran for his continuous interest and encouragement, and Mogens H0gh 
Jensen for interesting discussions of hydrodynamics. 
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